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SUMMARY AND ABSTRACT 


A subsonic flutter analysis capability has beer developed for NASTRAN, 
and a developmental version of the program has been installed on the CDC 6000 
series digital computers at the Langley Research Center. The flutter analysis 
is of the modal type, uses doublet lattice unsteady aerodynamic forces, and 
solves the flutter equations by using the k-method. Surface and one-dimensional 
spline functions are used to transform from the aerodynamic degrees of freedom 
to the structural degrees of freedom. Some prelininary applications of the 
method to a beamlike wing, a platelike wing, and a platelike wing with a folded 
tip are compared with existing experimental and analytical results. 


INTRODUCTION 


The available standard level of the NASA structural analysis computer 
program (NASTIAN) can be used to solve flutter problems by using the "direct 
input matrix" feature of the program to add the required unsteady aerodynamic 
force matrices to the appropriate structural matrices and solve the resulting 
eigenvalue problem. This procedure is inefficient and is not routinely used 
by aeroelasticians. However, since its first public release in 1970, NASTRAN 
has proven to be a very useful tool to many persons interested in flutter, but 
this use has been limited to using the program to calculate the structural 
modes euid frequencies that are required as input to separate special-purpose 
flutter analysis con^niter programs. This use of NASTRAN by aeroelasticians 
has created seme interest in incorporating a flutter analysis capability in 
NASTRAN. At the first NASTRAN Users' Experiences Colloquium (ref. 1) a paper 
(ref. 2) was presented that described the results of a design study for a 
ccanplete NASTRAN aeroelastic analysis capability. By using this design study 
as a guideline, the NASA has sponsored the development of a subsonic flutter 
analysis addition to NASTRAN. 


The purpose of this paper is to describe this new flutter analysis capa- 
bility euid present some results from preliminary applications of the program. 
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The technique developed is of the modal type, uses doublet lattice unsteady 
aerodynamic forces, uses one-dimensionaD. and surface spline functions to trans- 
form from aerodynamic degrees of freedom to structural degrees of freedom, and 
solves the flutter equations by using the k-method. The program is in what 
might be termed a developmental form, has only been installed on the CDC 6000 
series digital computers at the Langley Research Center, and is not available 
for general release to the public. Results from preliminary applications of 
the program to a beamlike wing, a platelike wing, and a platelike wing with 
a folded tip are compared with existing analytical and experimental results. 


OBJECT I VLS AND GUIDELINES 


The basic steps required in a flutter analysis are sho’m in the block 
diagram presented in figure 1. A characterization of each step is shown on 
the right in the figure. The overall objective of the NASTRAN subsonic flutter 
analysis was to provide a fully automated means for proceeding through these 
steps in an efficient manner to determine the flutter characteristics of 
complex structural and aerodynamic configurations. The development was con- 
strained to the use of existing, proven state-of-the-art techniques. There- 
fore, the major effort was to assemble the selected procedures into the NASTRAN 
environment. One of the most significant guidelines was that there would be 
no constraints imposed on the structural Idealization by aerodynamic considera- 
tions, and that the aerodynamic idealization would be made totally independent 
of structural modeling considerations. That is, the structure can be repre- 
sented by an optimum selection and arrangement of structural elements and 
degrees of freedom, and the aerodynamic characteristics can be determined by 
an optimum selection of aerodynamic degrees of freedom. This guideline 
dictated providing a very general capability for the structural- aerodynamic 
interface which is required to transform the aerodynamic degrees of freedom to 
the structural degrees of freedom. Additional guidelines were that the tech- 
nique should be easy to use and that the input data requirements associated 
with the unsteady aerodynamics and flutter solution and the format of the out- 
put results be in a form not totally unfamiliar to aeroelasticians. Another 
guideline that should be mentioned is that, vdiere practical, the new proce . .res 
required for flutter analysis would be made as general as possible so tha' the 
basic capability can be easily expanded, if so desired at a later date, to 
pccommodate additional aerodynamic theories, flutter solution procedures, and 
so forth. Naturally, it was required that the flutter analysis be compatible 
with the existing NASTRAN general structural capability and contain such 
existing features as the restart capability. Further, it wab required that the 
flutter auialysis be incorporated into a standard level version (level 15.1 was 
chosen) so that the NASTRAN program which contains the flutter analysis will 
also have all the other basic capabilities. 


METHOD IMPLEMENTED 


A modal flutter analysis method has been implemented in NASTRAN. The set 
of linear equations of motion that must be solved to determine the flutter 
condition may be expressed in matrix notation in the following form: 

508 
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M = generalized structural mass 
K = generalized structural stiffness 

Q = generalized unsteady aerodynamic force (function of Mq and k) 
b = reference length 
Mq = Mach nTJinber 
k = reduced frequency, btu/v 
V = velocity 

uji * generalized modal coordinate 
p = fluid density 
A = complex eigenvalue 

NASTRAN already contains the capability of generating the generalized 
|ass and stiffness matrices required by equation (l) but does not contain any 
Internal aerodynamic force capability. So one of the major tasks was to add 
le required unsteady aerodynamics. Since an important objective was to be 
*ble to analyze the most general aerodynamic configurations possible, the 
ioublet lattice unsteady aerodynamics method was selected for inclusion since 
[his method is applicable to a broad range of configurations. The flutter 
olution method implemented was the k-method which is the one most commonly 
sed in flutter analysis. A modal formulation was chosen for two reasons. The 
irst reason is that this is standard practice; the second reason is that the 
rder of the final matrix equations that must be solved is relatively small, 
le aerodynamic- structural interface is accomplished by the use of one- 
limensiciial and surface spline functions. 


The k- Method of Solution 

The k-method of flutter solution requires the repeated solution of equa- 
;lon (1). The aerodynamic forces are functions of the three parameters, 

i ensity, Mach number, and reduced frequency. To solve equation (l) values of 
wo of the parameters, usually density and Mach number, are held constant, and 
le eigenvalue equation is solved repeatedly for different values of reduced 

I requency. Ttie way equation (l) is developed, the damping, velocity, and 
requency of the system can be determined from the eigenvalues by using the 
elationships 
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Since the flutter point is on the boundary between stable (damped) and unstable 
(divergent) sinusoidal oscillations, the flutter condition occurs for the 
psirticular combination of parameters that causes the damping to equal zero 
(g = 0). The flutter velocity is usually determined by graphically plotting 
the damping versus velocity (g - V plots) obtained for each solution of the 
eigenvalue problon. A number of loci, equal to the order of the problem, will 
be obtained. The curve which crosses the g = 0 axis at the lowest value of 
velocity determines the critical flutter condition. The k-method implemented 
in BASTRAN includes the generation of both damping and frequency versus veloc- 
ity plots (f - V plots). Also, the capability is provided for selecting any 
one of the three aerodynamic parameters as the one to be varied. 


Unsteady Aerodynamic Theory 

The unsteadj- aerodynamic theory implemented in the NASTRAN flutter analy- 
sis is the subsonic doublet lattice method (ref. 5 )» Of the available proven 
theories, this technique is probably the most general in that it can be applied 
to multiple nonplanar mutually interfering lifting surfaces and can be used 
to calculate body-lifting surface interference effects. The doublet lattice 
method adapted for NASTRAN use is similar to that described in references h 
and 5 * The program described in these references includes slender-body aero- 
dynamics to calculate body, or fuselage, forces but this feature has not been 
included in NASTRAN although the work required to implement body forces has 
been determined. 

The doublet lattice method requires that the aerodynamic surfaces be sub- 
divided into a grid of trapezoidal boxes. An example box arrangement is 
illustrated in figure 2. The analyst is required to specify the box arrange- 
ment subject to certain geometric constraints. For example, two of these 
constraints are that the boxes must be arranged in streamwise coltunns parallel 
to the free stream and that surface discontinuities such as fold lines must 
lie on box boundaries. The geometric constraints on the box arrangement are 
not severe and provide sufficient latitude to model adequately very general 
configurations. For the unsteady flow case, a spanwise line of acceleration 
potential doublets is placed at the one-quarter- chord station of each box. The 
doublets are related to pressure and hence to the force on each box. An aero- 
dynamic influence coefficient matrix is generated \rtiich relates the force on 
the boxes to the downwash on the boxes. The force acts at the one-quarter- 
chord point and the downwash point is the three- quarter- chord point. Both of 
these points are at the box midspan station. Typical force and downwash 
points are shown in figure 2. The downwash is a function of the streamwise 
slope and the vertical displacement normal to the boxes. Each box may be 
thougiit of in the context of being a finite element with the degrees of freedom 
(deflection at one-quarter-chord point, and deflection and slope at three- 
quarter-chord point) defined at two different points within each box. In the 
NASTRAN flutter development, it was decided that it would be desirable to have 
only one aerodynamic grid point for each box. The point selected was the 
center of each box. A transformation is used to convert the force and down- 
wash at the one-quarter and three-quarter- chord points of each box to corre- 
sponding forces and downwashes at the centers of each box. Therefore, there 
is one aerodynamic grid point for each box. 
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Structural-Aerodynamic Interface (Geometry Interpolation) 

One of the most significant features of the NASTRAN flutter analysis is 
the geometry interpolation capability that provides for the interconnection of 
the aerodynamic and structural models of the system. Since a very general 
capability is provided for the structural -aerodynamic interface, the structural 
model can be that best suited from structural considerations alone, and the 
choice of aerodynamic model is dictated by aerodynamic considerations alone. 

The geometry interpolation provides a transformation from the aerodynamxc 
degrees of freedom to the structural degrees of freedom. This transformation 
is accomplished by the use of one- dimensional and surface spline furictions. 

(See refs. 6 and 7*) The traditional one- dimensional spline has been general- 
ized to include torsional rotations in addition to bending deformations. Since 
these functions are based on the small deflection equations of infinite beams 
and plates, respectively, they are very good for the interpolation of the 
deformations of general structural systems. If the structure is expected to 
behave like a beam as would be the case for a high-aspect-ratio jet transport 
wing, the one-dimensional spline would be used; if the structure is expected 
to behave like a plate, say a low-aspect-ratio wing, the surface spline would 
be the appropriate choice. The use of combinations of the two splines is 
permissible and would be applied, for example, to a complete aircraft where 
the fuselage had the character of a beam and the wing was expected to exhibit 
platelike behavior. 


Aerodynamic Force Interpolation 

The k-method type flutter solution requires the solution of the flutter 
eigenvalue problem many times so that a relatively closely spaced sequence of 
points can be determined to make the g - V plots since the behavior of the 
loci of roots on the plot can often be quite complex and lead to misinterp’^e- 
tation of the results. Since one of the most expensive parts of a flutter 
analysis is the determination of the unsteady aerodynamic forces, it is 
desirable to actually calculate the aerodynamic forces for a minimvun number of 
values of the Independent aerodynamic parameter, Mach number, or reduced fre- 
quency. Fortunately, experience has shown that although the behavior of the 
solutions of the flutter equations as displayed on a g - V diagram may be 
complex, the variation of the aerodynamic forces with reduced frequency or 
Mach number is generally smooth and well behaved. Consequently, it has become 
more or less standard practice in aeroelasticity to evaluate the aerodynamic 
forces at a relatively small number of values of the Independent variable and 
Interpolate to determine the forces at additional values of the Independent 
parameter. This interpolation is relatively Inexpensive when compared to the 
cost of actually calculating the aerodynamic forces and results in the loss of 
very little accuracy. Aerodynamic force interpolation has been included in the 
NASTBAN flutter analysis. Both one- dimensional and surface splines are used. 

If the flutter calculations are limited to a constant Mach number, the linear 
spline is used to Interpolate over a range of reduced frequencies. If a set 
of aeroc^amlc forces have bem determined at two or more Mach numbers, the 
surface spline is used to Interpolate to intervening Mach numbers. Experience 
with the one-dimensional spline has shown that it ^s very good for aerodynamic 


interpolation. However, there are some indications that the accuracy of the 
surface spline technique, although it is satisfactory, is not as good as the 
linear spline. This is probably caused by the fact that the character of the 
three-dimensional behavior of the aerodynamic forces is not p^.ateli'ke. 


FUJTTER ANALYSIS RIGID FORMAT 


The assembly of the components of the flutter analysis into a NASTRAN 
rigid format (labeled Rigid Format J+5) required the use of many existing 
functional modules, the modification to a few existing modules, and the devel- 
opment of six completely new modules. An annotated block diagram of the new 
rigid format is presented in figure 3. The structural analysis section is 
essentially identical to existing Rigid Format 10 (Modal Complex Eigenvalue 
Analysis) down to the point of ccanplex eigenvalue analysis. The existing 
module PLOT was modified to accommodate plotting of the aerodynamic geometry. 
Both undeformed and deformed plots are available. Changes were made to the 
XYTRAN and XYPLOT modules for the purpose of making g - V and f - V plots. 
An upper Hessenberg method of complex eigenvalue extraction was added to 
module CEAD since this procedure is better suited to the requirements of 
flutter analysis than the two methods already available. 

The completely new modules are the Aerodynamic Pool Distributor (AFD), 
Geometry Interpolation (gI), Aerodynamic Matrix Generator (AMG), Aerodynamic 
Matrix Processor (AMP), Flutter Analysis Phase 1 (FAl), and Flutter Analysis 
Phase 2 (FA2). Module APD forms tables of aerodynamic data, defines the 
boundaries of the aerodynamic elements, sind locates and orients displacement 
components at aerodynamic grid, or control, points. Module AMG evaluates the 
aerodynamic Influence coefficient matrix at specified values of Mach number 
and reduced frequency, and determines the transformations needed to convert 
these matrices from the points required by the doublet lattice theory (one- 
quarter and three-quarter box chord stations) to the center of the aerodynamic 
boxes. Module GI generates the transformations required to give the struc- 
tural displacements at the center of the aerodynamic boxes in terms of the 
deformations at the structural grid points. The AMP module calculates the 
generalized aerodynamic force matrices by using the mode shapes determined in 
the structural part of the rigid format (READ module), the aerodynamic matrices 
determined in AMG, and the transformation infoimation calculated in GI. The 
module FAl prq;>8kres the modal matrices for complex eigenvalue extraction by 
module CEAD. Also, the interpolation of the aerodynamic forces is carried out 
in this module, if a solution is required for a combination of parameters for 
which the generalized aero^mamic matrices were not determined previously. 

The module FA2 gathers data for reduction and presei^tation. For example, the 
velocity and frequency are determined from the eigenvalues calcu** ' ted by CEAD, 
and a line of printer output is prepared for each loop throu^ the flutter 
solution. The three modules FAl, CEAD, and FA2 are in a loop within the rigid 
format. This loop is repeated until solutions have been obtained for all the 
reduced frequencies, Mach numbers, and densities requested. 
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PBELIMINARY APPLICATIONS 


The NASTRAN flutter analysis has been applied to some simple geometric 
configurations. The results of three of these applications, a bearalike wing, 
a platelike wing, and a platelike wing with a folded tip, are presented in 
this section. The NASTRAN recults are compared with other available analytical 
results and experimental data. Some discussion of the features of the NASTRAN 
analysis is included with the discussion of the applications. 

The first application is the 15* swept wing shown in figure 4. Additional 
information concerning this wing may be found in references 8 and 9- This 
model was essentially a swept beam, and the NASTRAN structural model used 
consisted of 10 BAR elements as shown in the figure. The aerodynamic model 
consisted of 24 boxes arranged in six spanwise divisions of four chordwise 
boxes each. Unlike the requirements of the structural part of NASTRAN where 
the coordinates of each structural grid point are required input, a large 
number of aerodynamic boxes (and aerodynamic grid points which are located at 
the center of each box) are generated from a minimum amount of information. 

The aerodynamic boxes are assembled into panels, or groups, where each panel 
contains several boxes. For the beam example, all of the boxes belonged to a 
single panel. Only a single bulk data card (actually a parent card plus one 
continuation card) was required to define the aerodynamic boxes for this 
exaaple. For each group, the only information required is the coordinates of 
the inboard and outboard leading-edge corners of the panel (points marked "a” 
and "b” in the fig.), tie inboard and outboard chords (indicated by C]^ and 
C2 in the fig.), the number of chordwise boxes, and the niunber of spanwise 
boxes if the boxes are to be equally spaced. If the boxes are not to be 
equally spaced, then the desired spacing is provided in terms of fraction chord 
and spem divisions. The boxes for this example are equally spaced. Also, note 
that different coordinate systems were used to define the structural and aero- 
dynamic models. A one-dimensional spline function was used for interpolation 
in this example. Presented in figure 5 are the results of the NASTRAN calcula- 
tions for this wing at a Mach number of 0.45 and a density of I.I85 kg/m5. 

Three modes were used in the analysis. The results are presented in the form 
of a g - V plot where only the critical root is shown. The circle symbols 
indicate the calculated points. The calculated flutter speed is determined by 
the point at which the line faired through the symbols crosses the g « 0 
axis. Indicated on the figure, in addition to the NASTRAN result, are the 
experimental flutter result from refereitce 8 and the calculated flutter result 
from reference 9 which were obtained using linearized lifting- surface theory. 
The NASTRAN calculated velocity is in good agreement with the experimental 
value. The calculated flutter speed from reference 9 about 5 percent lower 
than the NASTRAN calculated value. The agreement with respect to flutter 
frequency is not so good. 

The wing geometry, structural model, smd aerodynamic model for the plate- 
like wing are presented in figure 6. Copies of NASTRAN computer-generated 
plots of the structural and aerodynasiic models are presented in figure 7. The 
structural model consisted of 36 quadrilateral plate elements (QUAD2); the 
aerodynamic model consisted of 50 boxes, 10 spanwise divisions of unequal 
spacing, and five equally spaced chordwise boxes. As was done for the beam 
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model, the entire v.' ; • made ui5 a single aerodynamic panel. The surface spline 
was used to perfr’tr .''.e required structural-aerodynamic interface for this 
example. The ce‘ ^*ed g - V curve is presented in figure 8, and only the 

critical locus 'nts is shown. These results are for s. Macn nianber of 0.80 

and a density - • .' 00 kg/cu m. Four modes were used in the o-naiysis. The 

solid syml.ri.' ;• ’;„*ate calculated values for which the generalized aerodynamic 
forces we" . ; . : .iated. The open symbols indicate results obtained by using 

interpola: ; -yireralized aerodynamic forces. The calculated and interpolated 
results to lie on the same curve and could not be distinguished from one 

another had not different symbols been used. Tabulated on the figure are tht, 
N/.STRAN ( aiculated flutter speed and frequt.icy, and some unpublished analytical 
results, .'ilso included in the table are some NASTRAN calculated results for 
an aerodynamic model that had eight equally spaced chordwise boxes and the 
same spanwise arrangement shown in figure 6 for the 50-box case. The unpub- 
lished analytical results were obtained by using a doublet lattice computer 
program similar., but not identical, to the one modified for NASTRAN use. The 
surface spline was also used for the structural-aerodynamic interface in 
obtaining the unpublished result. The experimentally determined model natural 
frequencies were used to determine the generalized stiffnesses used in the 
unpublished results. Since the measured frequencies did noc agree precisely 
with the calculated frequencies, some of the 7 -percent difference between the 
two results may be attributed to this frequency difference. However, the 
results are still in good agreement. The two NASTRAN calculations gave 
essentially the same results. 

The final application to be discussed is a platelike wing with a folded 
tip. A photograph of this model is presented in figure &nd the geometry, 
structural model, and aerodynamic model are presented in figure 10. The tip 
fin is inclined with respect to the wing by 6o*. Copies of NASTRAN generated 
counter plots of the structural and aerodynamic elements are presented in 
figure II. The wing portion of this model was the same as the platelike wing 
previously discussed, and this portion was modeled in the same fashion as the 
plate wing (36 QUAD2 structur.,1 elements and 50 aerodynamic bores comprising 
one aerodynamic panel). An cuiditiorjil 6o QUAD2 structural elements were used 
to model the folded tip. The folded tip was a separate aerodynamic panel and 
was cotDposed of a total of 50 boxes that were arranged into five equal chord- 
wise divisions and 10 unequal spanwise divisions as indicated in the figure. 

One provision provided by the program is that there may or may not be aero- 
dynamic interference, or coupling, between boxes located in different panels, 
or groups, dq>ending on the user to make the selection. This f-'ature allows 
for the cession of coupling when it is known to be unimportant and thereby 
reduces the time required to compute the aerodynamic matrices, or allows for 
the independent Investigation of aerodynamic interference effects. In the 
present exanqole, aerodynamic coupling between the wing panel and the tip panel 
was included. The surface spline option was used to perforo the required 
aerodynamic- structured interface. Four different spline functions were used, 
two for each aerodynamic panel. The i:.t.erpolation for the 25 inboard wing 
aeiodynasde boxes used one spline function, and the 25 outboard boxes used 
another spline function. The same type of arrangement was used for the tip 
fin. Since the analyst specifies the structural grid noints that are to be 
used for Interpolating for each aerodynamic box, it is not necessary that a 
single spline function be used for each aerodyn^c panel. 


The results of NASTRAN calculations for a Mach nuniber of O. 9 O and a 
density of O. 86 I kg/m^ are presented in figure 12 in the form of the g - V 
plot for the critical eigenvalue. Four modes were used in this analysis. The 
data obtained by using calculated generalized aerodynamic forces are indicated 
by the solid synbols in the figure, and the results using interpolated general- 
ized aerodynamic forces are indicated by the open symbols. The comparison of 
the results using calculated generalized aerodynamic forces with those obtained 
using interpolated forces indicates that they all lie on the same g - V curve. 
Also tabulated on the figure are an unpublished calculated result and an unpub- 
lished wind-tunnel experimental result. The unpublished calculated result was 
obtained in a fashion similar to that previously described for the platelike 
wing example. The two calculated results are in good agreement with respect 
to both flutter velocity and frequency. The experimental flutter velocity is 
about 9 percent lower than the hASTRAN calculated value. Both calculated 
flutter frequencies are somewhat higher than the experimental frequency. 

In discussing these three applications, some mention has been made of the 
simplicity of the input data requirements associated with the aerodynamics 
portion of the NASTRAN program. This point is somewhat dramatically indicated 
by the fact that for the wing with tip fin case, of a total of 401 bulk data 
cards used, only 28 were directly associated with the aerodynamics or flutter 
solution. 

Since the NASTRAN flutter analysis is relatively new, its efficiency has 
not been fully evaluated nor have all of its potential options been exercised. 
However, it is of interest to examine some of the central processing unit (CKJ) 
computer times required by some of the individual functional modules for u 
program execution. Presented in figure 13 is a listing of CPU times for the 
CDC *1600 computer obtained for the wing with the folded tip fin. In this case, 
five modes %rere calculated by the real eigenvalue module, and the four lowest 
modes were used in the flutter analysis. The generalized aerodynamic forces 
were determined at three values of reduced frequency and interpolated to two 
additional values so the flutter eigenvalue problem was solved five times. 
Additional information describing this example is shown on the figure. Also 
included on the figure are the total CHJ time, the peripheral processor time 
(CPU), and calls to the operating system (O/S calls). 


CONCUUDIRS REMARKS 


A subsonic flutter analysis capability has been developed for NASTRAN. 
This flutter analysis is of the modal type, uses doublet lattice unsteady 
aero^niamie forces, and solves the flutter equations by using the k-method. 
One- dimensional and surface spline functions are used to transform from aero- 
dynamic degrees of freedom to structural degrees of freedom. This capability 
has been incorporated into a version of NASTRAN, and this version has been 
installed on the CDC 6OOO series computers at the Langley Research Center. 
This version is in a developmental stage and is not now available for general 
release. In this paper, a general description of the new flutter analysis 
rigid format las been presented. Results of some preliminary applications of 
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the NASTRAN flutter analysis to a heamlike wing, a platelike wing, and a plate- 
like wing with a folded tip have been presented, and these results compared 
with existing experimental and analytical results. 
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Figure 2.— Sample arrangement of doublet lattice aerodynamic boxes. 
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STRUCTURAL ANALYSIS 
Determines generalized 
mass and stiffness, mode 
shapes, and frequencies 


EXISTING MODULES, 
READ, CHKPNT, ETC. 


Existing niodulos, arrv.ngement 
similar to standard le’/el 
Rigid F ormat IQ . 


Aerodynamic Pool Dis 
APD 


tr ibutor ^ 


New module, prep? res tables of 
aerodynamic drta and defines 
aerodynamic elements. 


Undeformed Plots 
PLOT 


Modified module, plots 
aerodynamic elements. 


FLUTTER ANALYSIS 
ADDITION 


6eome:ry Interpolator 
61 


Aerodynamic Matrix Generator 

AMG 


Aerodynamic Matrix Processor 
AMP 


Flutter Analysis Phase 1 
FAl 


New module, generates transformation 
from aerodynamic to structural 
degrees of freedom. 

New module, calculates a rodynamic 
influence coefficients for specified 
reduced frequencies and Mach numbers. 


New module, calculates generalized 
aerodynamic forces. 


New module, prepares modal matrices 
for complex eigenvalue extraction 
and performs aerodynamic interpolation 


Reduced frequency, 
Mach number, and 
density loop 


Complex Eigenvalue Analysis 
CEAD 


Solution Set Data Recovery 
VDR 


Modified module, upper Hessenberg 
solution routine added. 


Existing module, prepares eigenvectors 
for output. 


Flutter Analysis Phase 2 
FA2 

— — ~<rMore cases?^ 


New module, gathers flu^ ter results 
for reduction and presentation. 


XYTKAN, XYPLOT 


Modified modules, g-V and f-V plots. 


Data Recovery 
DDRl, SDRI, S0R2 


Existing modules, transforms results 
to physical coordinates, determines 
forces and stresses. 


Deformed Plots 
PLOT 


Modified module, plots flutter mode 
shapes. 


Figure 5 *— Block diagram of flutter analysis rigid format. 
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Figure 4.- Geometry, structural modeling, and aerodynamic modeling of beamlike vang. (All 

linear dimensions are in centimeters. ) 







Figure 6.— Geometry, structural modeling, and aerodynamic modeling of platelilre wing. (All 

linear dimensions are in centimeters. ) 
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Figure 9- — Photograph of platelike wing with folded tip fin 
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